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Abstract 

On the basis of the Friedberg-Lee model we formulate a semi- 
classical transport theory to describe the phase-space evolution of 
nucleon-nucleon collisions on the quark level. The time evolution is 
given by a Vlasov-equation for the quark phase-space distribution and 
a Klein-Gordon equation for the mean-field describing the nucleon as 
a soliton bag. The Vlasov equation is solved numerically using an ex- 
tended testparticle method. We test the confinement mechanism and 
mean-field effects in 1 + 1 dimensional simulations. 
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1 Introduction 



In this paper we extend the estabhshed transport theories, which describe 
dynamics of heavy-ion-coUisions on the basis of nucleons - 0), to a 
theory whose basic ingredients are quarks which are the relevant degrees of 
freedom at bombarding energies from some GeV up to some TeV. 

The fundamental theory for the interaction of quarks is QCD. There are 
derivations of transport equations for quarks and gluons based on QCD in 
the literature 0, |TD[- The problem is the complexity of the resulting 



equations which makes them rather tedious for practical applications. All 
other approaches start with some approximations to QCD. For bombarding 
energies larger than ^ 100 GeV QCD can be treated perturbatively. This 
is the motivation for many parton cascade models (Fritjof [|ll[], Venus [|l! 
RQMD 0, HIJING H, parton cascade [ll 



We are interested in the low energy regime of a few GeV where the 
QCD coupling constant is too large to treat QCD perturbatively. The non- 
perturbative QCD effects can be modeled by a so-called mean-field, which 
provides the confinement and governs the dynamics of the quarks. 

For the description of properties of the nucleon we have many well- 
established static quark models at our disposal. As a starting point for a 
dynamical theory we need a model, which generates dynamically a surface 
and is able to simulate absolute confinement. The last point is essential be- 
cause there will be excitations of quarks in a collision which could lead to 
deconfinement if the confining potential is finite. The simplest model which 
fulfills our requirements is the well-known Friedberg-Lee-soliton model with 
a field-dependent coupling constant (for a review see [0). Starting from 
this model we will derive transport equations for quarks moving in a mean- 
field. In this paper we study the resulting model in 1+1 dimensions to check 
the numerical methods and the behavior of the model in nucleon-nucleon 
collisions. 

Zhang and Wilets have derived transport equations based on the 
Nambu-Jona-Lasinio model in order to estimate chiral symmetry effects in 
heavy-ion collisions; this is conceptionally close to our work. However, be- 
sides the absence of confinement in their model, these authors do not actually 
perform dynamical simulations. 

This paper is organized as follows: In section ^ we review the static 
Friedberg-Lee model and present results in one spatial dimension. In section 
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1^ we derive equations of motion for the phase-space evolution of the quarks in 
this model. Section || discusses the extension of the usual testparticle ansatz 
to include particle and antiparticle degrees of freedom. The initialization 
of a stationary nucleon in the semiclassical approximation is discussed in 
section |^. First simulations of nucleon-nucleon collisions in 1+1 dimensions 
are presented in section ^ Finally, we summarize and conclude in section ^ 



2 The Friedberg-Lee Model 



The Friedberg-Lee soliton model in its basic version was formulated in 1977 
and 1978 by Friedberg and Lee []T^,|[T^. The quark-quark interaction in 
this model is mediated by a selfinteracting scalar a field. The scalar field 
is interpreted as a summation of all nonperturbative gluonic interactions 
between the quarks. Because of its color neutral nature, the cx-field can only 
model many-gluon-exchange. Therefore, the model has been extended to 
include absolute color confinement by introducing a color-dielectric function 



[ p7[l . In this paper we use the following simplified version of this model p8 



C = i^ (7^,9^ - mo) vl> - + - {d,af - U{a), 



(1) 



which describes quarks with a rest mass mo coupled to a scalar field a. This 
coupling together with the mass term leads to an explicit breaking of chiral 
symmetry. For the effective coupling to the scalar field we use the form given 
by Fai et al. [EH] (a similar coupling is discussed in 
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It is shown in that this form of the effective coupling gives a good ap- 
proximation to the effects of the gluon field in the chirally invariant color 
dielectric extension of the soliton-bag model. It guarantees absolute con- 
finement, because the quarks acquire an infinite effective mass mo + geff{(j) 
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in the vacuum. The 9- function in (|^) guarantees that the effective quark 
mass is larger than the quark rest mass and always positive (^f > 0), which 
is essential for our semiclassical treatment of this model. Therefore we start 
with the Lagrangian (|l]) to formulate a transport theory. 

For the scalar field cr, a nonlinear self-interaction U (a) is assumed, which 
is necessary to allow for solitonic solutions of the field. This nonlinear po- 
tential is parametrized as 

= f + + + • (5) 

The free parameters of the model Qq, a, b, c, B can be adjusted to reproduce 
the basic properties of the nucleon, namely mass, RMS radius, magnetic 
moment and the ratio Qv/ Qa- 

Solutions of the Friedberg-Lee model have been extensively studied dur- 
ing the last decade [|l6l , ||2^ - [|27 . 



In this paper we restrict our studies to 1 -|- 1 dimensions. In order to 
ascertain that the essential properties of the model are still present in 1 -(- 1 
dimensions we have compared the solutions of the Friedberg-Lee model in 
1 spatial dimension with the 3-dimensional results from the literature. It 
has turned out, that the results are essentially identical. In order to show 
this, we present a 1-dimensional result for the model (|I|) with the effective 
coupling given by @. In one space-dimension the ansatz for the spinless 
quark spinors is 

(«) 

which gives the following equations of motion in the usual mean-field ap- 
proximation 

— = -{e + ga + mo)v, (7) 



dx 

dv 
dx 



ga - mo) u, (8) 



= N'4^iu^-v^) , (9) 
dx^ da da 
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with the occupation number = 3 for the valence orbitaL A typical solution 
for the parameters a = 222.91, b = —5347.05, c = 38610, go = ^ with a total 
energy of 931MeV and an RMS radius of 0.56/m is shown in figs. ^ and 0. 



3 Derivation of the Transport Equation 

The aim of this section is to derive equations of motion which describe the 
phase-space evolution of a moving nucleon. We therefore follow the spirit of 
the well known RBUU model 0, which has been applied very successfully 
to heavy ion collisions 0. The following steps are similar to the derivation 
in pOl ^. Starting with a static model for the nucleon, the Friedberg-Lee 
model, we derive an equation for the time evolution of the Wigner-function, 
which is the quantum mechanical analog of the classical phase-space density. 
The Wigner-function is defined as 

W'^-'^{x,p) = / d^R e-'P"«'^^'(x + ^)® ^'■^■(x - ^) , (10) 

where the indices rj, label internal degrees of freedom, such as color, flavor, 
etc; in the following they are suppressed. Densities and particle number can 
be easily calculated by integrating over the Wigner-function : 



p{x) = J d^ptr{^oWix,p)) (11) 
p{p) = J d^xtr{^oW{x,p)) (12) 
N = J d^x J d^ptr{-foW{x,p)) . (13) 

In general, expect at ion- values of one-particle operators O are given by 

(^O) = Jd^xJ d^p tr{dW{x,p)) . (14) 
The equation of motion for the Wigner-function can be derived by calculating 

[7^ {d>^ - 2tp'^)] W{x,p) = -^Jd'R e-''P^^''^{x,)®^,d:^^{x2) , (15) 

where ^ ^ 

Xi = X + — , X2 = X - —. (16) 
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On the r.h.s. of equation (|1^) one can use the Dirac-equation of the Friedberg- 
Lee-model to replace 71,9^2^(^2) which gives after some algebraic transfor- 
mations: 

_2i e^M^.d^. [rrio + ga{x)] W{x,p) . (17) 

The derivative (9^ acts on the Wigner-function and d!^ acts on the field a{x). 
The factors h are inserted explicitly. Up to now, equation (0) together with 
the equation for the mean field 

d.d^a _^ = ^-lfd.p trWix, (18) 
da da J 

is equivalent to the equations of motion for the Friedberg-Lee model in the 
mean-field approximation. In the semiclassical approximation, one expands 
the Wigner-function and the exponential function in equation { [T^ ) in orders 
of h: 

W = Wo + inWi + --- (19) 



2 



This expansion is expected to converge the better the less the fields vary over 
the Compton wavelength of the quarks. In lowest order h this leads to the 
equation 

{^^p'^-m*)Woix,p)=0 , (20) 

with m* = mo + geff{cr{x)) 

This is the well known mass-shell constraint. In first order h one has: 

7m^^ + dPd^m*{x)] Wo{x,p) = 2 [7^^^^ - m*] W,ix,p) . (21) 

We now take the trace on both sides of this equation. If one requires baryon 
current conservation in lowest order in h, Wi has to fulfill a constraint so 
that the trace of the r.h.s. of equation (|2TD vanishes. We are then left with 
an equation of motion for Wq, the well known Vlasov-equation. 

All equations are so far derived for 3 + 1 dimensions. From now on we 
study the model in 1 + 1 dimensions in order to test its main features and the 
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numerical methods for its solution. The generalization to 3 + 1 dimensions 
will be discussed in a future work. In 1 + 1 dimensions we introduce Dirac- 
Matrices F^, which have to fulfill the usual anti-commutation relations: 

{T^,T,} = g^, /i,z/e{0,l} (22) 

with g^^= (^^ - 1 ) ' ^^^^ 
One possible representation for the F-matrices is 

One also defines the product 

r2 = FoFi=(^J J)' (25) 

which is analogous to 75 in 3+1 dimensions. 
The relations 

{F2,F4 = , Tl = I , Tl = -I , [/,andF^]=0, (26) 

hold for these matrices, where / is the 2x2 unit matrix. 

In the 1-dimensional notation the Vlasov-equation can be written in the 
following form: 

tr {[Todt + TA, + I{dtm*dE - d,,m*dpi)] Wo{x,p)) = 0, (27) 

with X = (t, x^) and p = {E,p^). 

This equation incorporates energy- and particle-number conservation; it 
approximates the quantum-mechanical solution only if the fields vary slowly, 
a condition which we will discuss later. 

In the following we deal with Wq only and suppress the subscript. In view 
of the mass-shell constraint eq. (|20|), one takes the following ansatz for the 
Wigner matrix: 

W{x,p) = {T^p^ + m*)f{x,p) , (28) 
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where f{x, p) is a scalar function of x = {t, x^) and p = {E,p^). This is equiv- 
alent to the usual spinor decomposition neglecting the pseudoscalar part. 
Insertion of this ansatz in equation ( jn\ ) finally gives the Vlasov-equation: 

'p^d^ + m*d^m*dP] fix, p) = 0. (29) 

The mass-shell constraint 

[T^p^ -m*)W{x,p) = (rX-m*)(rX + m*)/(x,p) = 

= (p^ -m*^) f{x,p) = (30) 



will be fulfilled by a using a Dirac 5-function 6 [p^ — m*^^ in f[x,p). For the 
numerical treatment of this equation we use the so called testparticle ansatz 
for the scalar phase-space distribution function f{x,p) jSlj : 



fix, p) = 6{p'~ m*') e{E) Y: Six - xM)6{p' - pl{t)) . (31) 

n 

The ^-function expresses the restriction to positive energy states. This ap- 
proximation will be discussed in detail in section |. Inserting this ansatz into 
the Vlasov-equation ( P5| ) shows that the so called testparticles have to move 
like classical particles according to Hamilton equations of motion: 

Xn = ^ (32) 

m* 

pI = -T^d.K ■ (33) 

with 

En = ^pi^ + m^^Xn) (34) 

4 Extended testparticle ansatz 

The usual testparticle ansatz ( |3TD which is restricted to positive energy so- 
lutions implies: 

dEtrW{x,p) = W,,alar{x,p^) = — -r^f{x,p^) (35) 

J dEtr{roW{x,p)) = Wbaryon{x,P^) = f{x,p^) . (36) 
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Therefore, at an initial time instant a given baryon-density can be repro- 
duced by using eq. (^61). The scalar density is then given by eq. (|35|) . In 
nuclear transport theories the rest mass of the testparticles is of the order of 
IGeV. In this case this ansatz works quite well. On the quark level the rest 
mass of the quarks is of the order of 5 — lOMeV, which is reflected in the 
fact that the lower component of the quantum mechanical wave functions is 
nearly as large as the upper component (see fig. A projection on nega- 
tive energy states is not negligible. It is therefore necessary to extend the 
testparticle ansatz in order to fit scalar and baryon density simultaneously. 
The importance of the scalar density is obviously clear because it determines 
the dynamics of the system through its coupling to the cr-field, whereas the 
baryon density determines the charge radius of the nucleon. 
Our suggested extension of (|3T|) is given by 

fix,p) = fP'''^%x,p)9iE) + (l - CJ;{x,p)) e{~E\ (37) 

where /poJ**'^'*^ represents the occupation of positive energy particle states and 
(1 — fnfg) the occupation of negative energy hole states relative to the filled 
Dirac-sea. With this ansatz the Wigner-f unction is written as 



2\E\ 



+ (rX + m*) 0{E)'^-^^^^^^ (l - d^t, x,p)) e{-E), 



with oj = sjp^ + m*^, and x and p denote the second component of the 2 
vectors x^^p^. The energy dependence of the on-shell distribution functions 
fit, x,p) is implicit. 

Performing the energy integration we end up with: 

+ -'°"-^"'^"'' (l-/„y«.x.ri) .(38) 
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where uj is the on-shell energy. Defining 

f^::''^%t,x,p) = f{t,x,p) (39) 

= f{t,x,-p) (40) 

this can be rewritten as 

fdE Tquj - Tip + m* Tquj + Tip - m* - 
y 2^W^(^,2;,p,£;) = 2;-^ f{t,x,p) — f{t,x,-p). 

(41) 

The meaning of / and / becomes more clear if one calculates the scalar 
and the baryon density: 

dEdp trW = j dp— (/ + /) 
Pb{x) = JdEdptr{ToW)^ Jdp{f-f). (42) 

While / und / add up in calculating the scalar density, their difference yields 
the baryon density. In view of this / can be interpreted as particle distibution 
function whereas / is the anti particle distribution function. 

For a given Wigner-function W f and / will then be constructed accord- 
ing to 

f{x,P) = ^ (^-^Wscalar + Wbaryon^ 

f{x,p) = ^ (^-^Wscalar - Wbaryon^ , (43) 



with 
W, 



'scalar ^ j trW{t, X, E,p) , Wbaryon ^ j dE tv {ToWit, X, E,p)) 

(44) 



and id = yp'^ + m* . 

For the numerical realization we initialize A^^i testparticles according to 
the distribution f{x,p) and N2 testparticles according to f{x,p), where 



with Ii = dxdp f{x,p) 

and . (46, 



All testparticles follow the Hamiltonian equations of motion (p2|), {poj}. This 
can be easily understood, because in a scalar field particles and antiparticles 
feel the same forces. The baryon density and the scalar density are given by 



1 ^ 



^ n=l 



1 A m* 



Psi^) = ^E^^i^-^n), (48) 
n=l 



where N = Ni + N2, N = Ni - N2 
and 6< 



1 if testparticle from / 
1 if testparticle from /. 



On a spatial grid with spacing Ax the 5-functions are evaluated as 

[ otherwise. 

With this extended testparticle ansatz, it is possible to reproduce a given 
baryon density and a scalar density independent of each other. In the se- 
miclassical approximation the phase-space evolution of the Friedberg-Lee 
model is determined by the equations of motion for the testparticles (^) 
and (^) and by the mean-field equation by 



5 A semiclassical nucleon 

We first tried to initialize testparticles as described in the last section to 
reproduce the baryon and scalar phase-space distributions calculated from 
the quantum mechanical solution (fig. ^). Performing the time evolution of 
the testparticles and the a-field, it turned out that the nucleon initialized in 
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this way is not a stable solution of the semiclassical equation of motion (|27|). 
This failure indicates that the conditions of weakly varying fields necessary 
for the derivation of the transport equation are actually not met. Similar 
difficulties show up when we calculate the semiclassical energy-distribution 
from the quantummechanical solution. The latter, is of course, a 5-function 
whereas the former turns out to be a broad distribution. 

We therefore drop now the connection between the semiclassical and the 
quantummechanical solution and from now on work consistently in the trans- 
porttheoretical framework. In doing so we take the potential function U{a) as 
an effective classical potential that incorporates all higher order corrections 
to the quantummechanical one and determine it by adjusting the properties 
of the semiclassical solution to empirical properties. 

First, we look for solutions of the stationary Vlasov-equation. A station- 
ary solution can be characterized by 

dtfix,p) = and dtui* = dt geffia{t)) = . (51) 

The Vlasov-equation reduces to 

E E 

— -dj{x,p) = —dpf{x,p) , (52) 

where E = \Jp'^ + m*^. Every function f{E) fulfills this equation. This 
means that every distribution function which does not explicitly depend on 
position and momentum, but only on the combined quantity energy is a 
stationary solution of the Vlasov- -equation. For the initialization of a stable 
nucleon we have to find a distribution which reproduces the properties of the 
nucleon. 

We start with a functional ansatz: 

ftotai{x,p) = f{E) + f{E) (53) 

where 

f{E) = /oe-'^^^ 

m = foe-'^' (54) 

are two Gaussians in energy with four free parameters; the amplitudes and 
the widths. The soliton solution is then calculated by starting with a Wood- 
Saxon shape for the a-field. With the Gaussian ansatz we can also calculate 
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energy and density contributions: 



Ps = / cipf (/o e--^ + /o e--^ 
PB = J dp {fo e-'^' - fo e-'' 

Equark = J dp J dxE (/q S'''^ 



1 (da 

2 [d^ 



E^ = Jdx \^U{a) + 

Etotal = '^Equark + -E'er, (55) 




with E = \Jp^ + m*^ 

The (T-field is then calculated selfconsistently with the scalar density. The 
parameters are varied to give reasonable values for the energy and RMS- 
radius of the nucleon. In principle also the parameters for the potential 
could be varied to find the optimal parameter set for this 1+1 dimensional 
case, but this was not done here. 

In fig. 1^ we show a typical solution. Baryon density and the cr-field are 
similar to the quantum mechanical solution. The baryon density is volume 
centered. The shape of the scalar density is quite different from that of the 
baryon density and surface-peaked, reminiscent of the SLAC bag solutions to 
the Friedberg-Lee soliton bag model |TB|. The surface-peaking is here a clear 



consequence of the strong increase of the coupling constant g towards the 
surface. The maxima at the surface are stabilizing, because the testparticle 
mass increases at the surface which results in a deceleration, meaning that 
the testparticles stay longer in the surface region. Fig. ^ and ^ show the 
stability of the semiclassical nucleon in the time evolution. The stability is 
almost perfect. 



6 The boost 

After the discussion of a stable initialization for a nucleon we now describe 
the Lorentz- boost of a given testparticle distribution, necessary for the 
preparation of the initial state of a collision. The Lorentz-boost has to 
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give the correct transformation properties for the baryon density, the scalar 
density, the cr-field and the Lorentz scalar distribution functions / and /. 

We initialize the field and distribution functions at a time t' = in the 
moving frame. The transformation is therefore given by: 

t' = 0, x' = - (56) 

7 

p' = ^{p-PE), E' = ^iE-(3p) (57) 

To calculate the cr-field and the distribution functions / and / at the time 
t' = 0, we need to know the configuration in the rest frame for different 
times t. In this frame, however, the nucleon is stationary, i. e. f{t) = /(O) 
and a{t) = a{0) for all times t. Therefore, we can set t = 0. 

The assumption that the distribution function is a scalar implies that 

f{x',p',t' = 0) = f{x,p,t = 0) = / (7x',7(p' + PE'),t = 0) (58) 

Inserting the testparticle ansatz: 

fix, p,0)=Y. ^(^ - a;n(0))5(p - p„(0)) (59) 

n 

we get 

f'ix\p', 0) = ^ 5 (7x' - x„(0)) 6 (7(p' + f3E') - p„(0)) . (60) 

n 

We express this sum using the coordinates in the moving system: 

fix', p', 0) = ^(1 - (3v^)6ix' - x^'iOMp' - p„'(0)) (61) 

n 

X 

with Xn = — and Pn = liPn - l3En) . (62) 
7 

In addition to the usual transformation of the testparticle coordinates and 
momenta often used in the literature we obtain a factor (1 — /3f„) in the 
testparticle sum which results from the transformation of the arguments in 
the 5-f unctions. 

As a check for the expression for the boost (|6lD , we calculate the trans- 
formed densities: 

ps'ix') = ( fix',p')dp' = Y.5ix-xM'^-[3vn)=lipBix)-P3iix)) 

n 

/m*' . m* 

—fix',p')dp' = Y.5ix-x^)^ = p,ix) . (63) 
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Here, ji denotes the baryon current, which vanishes for a stable nucleon at 
rest. The densities have the correct transformation properties. The antipar- 
ticle distribution f'{x',p') will be correspondingly transformed. 

If one inserts the testparticle sum ( |6lD into the Vlasov-equation, the 
resulting testparticle equations contain terms generated by the momentum 
derivatives acting on the factor (1 — To avoid this change, we use a 

new testparticle distribution with coordinates x,p, which have to fulfill: 

- = ^(1 - - a;„')<^(p' . (64) 

n n 

This can be guaranteed, if the new testparticles located at x,p are distributed 
according to the modified distribution function in the rest frame 

f{x,P) = {l-(3v)f{x,p) . (65) 

Finally, a moving nucleon has to be initialized in the following way: The test- 
particle coordinates are distributed according to ( |65D and the transformation 
to the moving coordinates is given by (0). The a-field remains unchanged 
with the transformed x coordinates. 

7 Collisions 

To simulate collisions of two nucleons, we initialize semiclassical nucleons 
as described in section |^ and transform them to the center of mass system 
by applying the boost of section ^. For each of the nucleons the center 
of the boosted phase-space distribution in momentum is shifted to pc = 
7/9mo; due to the small quark rest mass, it stays close to zero momentum. 
This means that the boosted phase-space distributions always overlap in 
momentum space around p ^ 0. The boosted phase-space distribution is 
symmetric in x direction. In momentum direction one half of the distribution 
is contracted and the other is stretched, because of the factor (1 — /3f„) in 
(pT]). This leads to a very asymmetric shape in momentum space. 

All simulations which are shown in this paper are calculated using the 
testparticle method. The testparticle propagation and the mean field time 
evolution are calculated using predictor-corrector techniques. As a check 
we also integrated the Vlasov equation numerically on a phase-space grid 
using an alternating direction implicit procedure (ADIP) which leads to the 



14 



same results. The direct numerical integration can be performed only in the 
1+1 dimensional case for reasons of computing time. For the extension to 
3 spatial dimensions the testparticle method is clearly superior for practical 
applications. All calculations are done in the c. m. system. In the simulations 
presented here we use 20000 testparticles for each nucleon. This relatively 
large number in the 1-dimensional case is necessary to obtain good statistics 
for the calculation of the scalar density. Due to the increasing coupling- 
constant at the surface of the soliton, the cr-field is very sensitive to small 
changes of the density in this region which requires good spatial resolution 
that cannot be obtained by the rather coarse-grained smearing methods. 

First, we have tested the stability of one moving nucleon. It is per- 
fectly stable for times (~ 20/m/c) larger than the typical collision times 

1/m/c). The time evolution of the quark distibution and the sigma field 
is shown in fig. ^ and fig. ^ respectively. 

We have calculated collisions for bombarding energies of 1 — 25 GeV in 
the laboratory frame. In the simulations we initialize two nucleons which are 
well separated in space; the propagation of a few fm until the collision takes 
place ensures that we have stable incoming nucleons. For a boost velocity of 
(3 > 0.9 the collision time is of the order of 1 fm/c. To be sure to observe 
stable outgoing nucleons, we continue the calculation for 5 — 10 fm/c after 
the collision. 

For very low bombarding energies {Tiab ~ IGeV) we observe a fusion of 
the two incoming bags. This is due to the slow motion of the nucleons which 
leads to a long overlap time. The testparticles are decelerated and loose 
kinetic energy which results in a gain of potential mean-field energy. The 
result is one extended bag with constant density inside. Due to the mixing of 
testparticles this bag does not break up into separate bags. The fusioned bag 
reaches a maximal extension and shrinks again converting potential energy 
back to kinetic testparticle energy. This behaviour is shown in fig. |^. The 
oscillation is undamped because there is no mechanism for energy loss in the 



model. We note that Schuh et. al. 32] estimate the potential for nucleon- 



nucleon scattering in the Friedberg-Lee soliton bag model and also end up 
with an attractive potential which may be responsible for the fusion that 
we observe for slowly moving bags. However, the surface energy in our 1+1 
dimensional case is different from the real 3+1 dimensional world so that it 
is not sure that the fusion effect will survive in 3 spatial dimensions. 

In collisions with bombarding energies Eiab from 5 — 25 GeV the nucleons 
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are transparent on the mean-field level. The phase-space evolution of the 
baryon distribution (/ — /) is shown in fig. |^ for a bombarding energy Eiab = 
15 GeV. Since the scalar field acts in the same way on testparticles and 
antitestparticles, the distribution functions / and / show the same behaviour. 

In fig. ^we display the corresponding a-field and the scalar source density. 
In the moment the two nucleons collide, the scalar density pushes the cx- 
field barrier down. The a-field overshoots to negative values (timestep t = 
1.1 fm/c), but due to the 6'-f unction in the effective quark-cr couphng this 
does not lead to negative effective masses. The field decouples from the 
density in the region where it is negative and oscillates up again (timestep 
t = 1.2 fm/c). In the following time-step {t = 1.3 fm/c) the cr-field comes 
actually very close to the vacuum value. Here the coupling to the quarks then 
becomes large as shown by the large values of the scalar density. The quarks 
thus become very massive and stay confined. Remembering the physical 
motivation for the ingredients of our model we may say that at this point the 
explicit gluon-exchange forces take over from the color background field. The 
correlation then swings back to the background field thus reestablishing the 
soliton. The velocity of the nucleons is too large to allow for any significant 
one-body dissipation during the short overlap time. The result is a dip 
in the density in the moment the nucleons have passed through each other 
(timestep t = 1.6 fm/c). The cr-field then comes close to the vacuum value 
and the two solitons separate (timestep t = 1.9 fm/c). If we continue the 
simulation, the small deformations, which can be seen in phase-space at 
t = 1.9 fm/c, vanish and we end up with two stable nucleons. The time 
evolution of the energy is shown in fig. ^. There is no significant exchange of 
energy from the testparticles to the mean-field. This is most probably due 
to the restricted geometry in the present version of our model. A realistic 
surface in a 3 + 1 dimensional calculation would clearly lead to a larger 
coupling between mean-field modes and the quark motion. 

We conclude from these simulations that the model is able to describe 
moving quark systems and absolute confinement in a collision of nucleons. 
Our model is in this sense a basis for extensions in the direction of direct 
quark-quark collisions, which we will put on top of the stable colliding nu- 
cleons we describe so far. 
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8 Summary and Conclusion 



In this paper we have formulated a transport theory on the quark level which 
includes nonperturbative aspects of QCD. This is done by including a soliton 
mean-field which governs the dynamics. As a starting point we chose a 
phenomenological quark model for the nucleon, namely the Friedberg-Lee 
model. We have derived equations of motion for the time evolution of the 
Wigner-function which leads in the semiclassical expansion to the well known 
Vlasov-equation combined with a mass-shell constraint. For the Wigner- 
function we have used a testparticle ansatz which results in classical equations 
of motion for the testparticles. 

Considering different extensions to the original Friedberg-Lee model it 
has turned out that the formulation with an effective quark-cr coupling, which 
comprises non~color-singlet many-gluon-exchange effects, can be used to 
model absolute confinement on the testparticle level. Inside the dynamically 
generated bag the quarks are nearly massless ('asymptotic freedom'). 

The usual testparticle ansatz was extended including also negative energy 
states which was necessary because of the small rest mass of the quarks. 
With the extended ansatz it is possible to describe scalar and baryon density 
independently of each other. 

For the initialization of a stable nucleon we have used a consistent so- 
lution of the stationary Vlasov-equation. It was shown that any distribu- 
tion function which depends only on energy and not explicitly on position 
and momentum is a solution of the stationary Vlasov- equation. After the 
discussion of the semiclassical nucleon we have described how to boost the 
initialized nucleon to give boosted distribution functions with the correct 
Lorentz transformation properties. 

Having all these ingredients, we have performed some first collisions in 
1-1-1 dimensions in order to test the numerical methods and the behavior of 
the model. It has turned out that on the pure mean field level the nucleons are 
totally transparent for high bombarding energies. The confinement is realized 
for each nucleon separately. For very low bombarding energies (< iGeV in 
the lab. frame) we observe a fusion of the bags with a following oscillation of 
the six quark bag. This may be an artifact of the 1-dimcnsional treatment. 

At the moment we are extending our model to 3+1 dimensions. This 
is essential to allow for non-central collisions and sidewards fiow. The next 
step is then to include a direct quark- quark collision term and cross sections 
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for particle production. 
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List of figures 



Figure 1: Baryon density (dashed), scalar density (dotted) and cr-field 
(solid) for a typical solution of the Friedberg-Lee model with effective cou- 
pling in 1 dimension. (The cr-field is in arbitrary units.) 



Figure 2: Upper component u (solid) and lower component v (dashed) of 
the quark spinor for the 1-dimensional solution of the Friedberg-Lee model 
with effective coupling 



Figure 3: Typical solution of the stationary Vlasov-equation with the ansatz 
described in the text: baryon density (dashed), scalar density (dotted) and 
(T-field (solid). The parameters arc a = 245.75, b = —5614.398, c = 38610, 
/o = 7.527, k = 2.0, /o = 0.7528, k = 5.0. The total energy is SAlMeV and 
the RMS radius 0.437/m. 



21 



Figure 4: Time evolution of the quark distribution of a nucleon at rest in 
phase-space. The initiahzation is described in the text, (time in fm/c). 



Figure 5: Time evolution of the scalar density (dashed) and cr-field (solid) 
for a nucleon at rest, (time in fm/c) 



Figure 6: Phase-space evolution of the quark distribution in a collision of 
two nucleons for a bombarding energy of IhGeV in the lab. system. The 
time is given in fm/c. 



Figure 7: Time evolution of the scalar density and the cr-field for a collision 
of two nucleons with a bombarding energy of IhGeV in the lab. system. The 
time is given in fm/c. 



Figure 8: Time evolution of the energy for a collision with bombarding 
energy of IGeV in the lab. system. Total energy (sohd), testparticle en- 
ergy (dashed), cr-field kinetic energy (dotted) and cr-field potential energy 
(dashed-dotted) . 



Figure 9: Time evolution of the energy for a coUisioon with bombarding 

energy of IbGeV in the lab. system. Total energy (solid), testparticle en- 
ergy (dashed), a-field kinetic energy (dotted) and cr-field potential energy 
(dashed-dotted) . 



22 



